Geometric and impurity effects on quantum rings in magnetic fields 
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We investigate the effects of impurities and changing ring geometry on the energetics of quantum 
rings under different magnetic field strengths. We show that as the magnetic field and/or the 
electron number are/is increased, both the quasiperiodic Aharonov-Bohm oscillations and various 
magnetic phases become insensitive to whether the ring is circular or square in shape. This is in 
qualitative agreement with experiments. However, we also find that the Aharonov-Bohm oscillation 
can be greatly phase-shifted by only a few impurities and can be completely obliterated by a high 
level of impurity density. In the many-electron calculations we use a recently developed fourth-order 
imaginary time projection algorithm that can exactly compute the density matrix of a free-electron 
in a uniform magnetic field. 
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I. INTRODUCTION 

The possibilities of utilizing the Aharonov-Bohm (AB) 
effect in future nanotechnology has stimulated signifi- 
cant recent progress in fabricating nanoscopic quantum 
ringsiS^ Since the capabilities to control the electron 
number in the ring and to modify its geometry are essen- 
tial for observing and exploiting new phenomena, sys- 
tems with these capabilities have attracted much theo- 
retical interest. 

Lorke and co-workers^ have applied self-assembly tech- 
niques to create InGaAs and GaAlAs/GaAs rings con- 
taining only a few electrons. The weak electron-electron 
interaction in these rings makes them most suitable for 
optical experiments, and the observed state transitions 
can be well explained with the single-electron spectrum 
of a parabolic ringji On the other hand, Keyser and co- 
workers^ have reached the strongly correlated regime by 
omitting the screening gate on top of a few-electron quan- 
tum ring fabricated from a GaAs/AlGaAs heterostruc- 
ture. They were able to observe fractional AB oscilla- 
tions with a period of $o/N, where $0 = h/e is the flux 
quantum and N is the electron number. Exact diagonal- 
ization of a few-electron Hamiltonian^ has shown that 
the electron-electron interactions break the degeneracy 
between the singlet and triplet states, leading to frac- 
tional oscillations. Similar results have been obtained 
within the Heisenberg modelS. and also by recent Monte 
Carlo calculations. 7 - Both have clarified the role of the 
electron localization in the fractional AB effect. The role 
of the (two-dimensional) width of the ring is, however, 
still unclear in the strong- interaction limit. 

From magnetotransport experiments in the Coulomb 
blockade regime one can infer the discrete energy lev- 
els of a quantum ring^ Moreover, these measurements 
have been performed for both circular and square ring- 
geometries; the latter corresponds to a chaotic Sinai bil- 
liard^ i.e., a circular antidot at the center of a square 
quantum dot. In such a symmetry-broken quantum ring, 
the increasing magnetic field induces regularity in the 



amplitude and position of the Coulomb peaks. These 
dots have been estimated to contain hundreds of elec- 
trons. The interactions were screened by a top gate, 
hence the simple single-electron picture provides a suf- 
ficient description of the spectrum^ 

We examine in this paper the energetics of circu- 
lar, square-shaped, and impurity-doped two-dimensional 
(2D) quantum rings containing up to N ~ 20 strongly 
interacting electrons. We shall focus on the measur- 
able quantities such as the chemical potentials, addi- 
tion energies, and the magnetization. We will iden- 
tify the quasiperiodic AB oscillations as well as differ- 
ent magnetic phases, and find, in agreement with the 
experiments, 2 that these become very similar between 
circular and square rings at large N and in high mag- 
netic fields. In these phases we find similarities to in- 
teger and fractional quantum Hall states of quantum 
dotsAiS We also carry out a statistical analysis of the 
addition-energies for quantum rings containing randomly 
distributed Coulombic impurities. Increasing the num- 
ber of impurities leads to a systematic phase shift in the 
AB oscillations. As a result of the electron-electron in- 
teractions, the high-disorder limit is characterized by a 
Gaussian-like addition-energy distribution. This shape 
indicates the disappearance of the AB oscillations. 

In this work, we solve the many-electron problem in a 
strong, uniform magnetic field by use of the spin-density- 
functional theory (SDFT). To solve the Kohn-Sham (KS) 
equations for thousands of impurity configurations, we 
use our recently developed fourth-order projection algo- 
rithm^ to determine the occupied KS orbitals. This 
algorithm is highly efficient, since the number of Fast 
Fourier Transforms used for solving the KS spectrum 
remains the same even in the presence of a magnetic 
field. This is because at its core, our algorithm is capable 
of exactly computing the density matrix of a free elec- 
tron in an arbitrarily strong magnetic field. Thus unlike 
other methods of solving the Schrodinger equation, the 
magnetic field part of the physics is hardwired into our 
algorithm. Also, instead of the usual slow-convergent, 




charge-mixing iterations, we update the charge densi- 
ties by our linear-response algorithm 12 with accelerated 
Newton-Raphson convergence. These significant algo- 
rithmic advances are not based incremental improvement 
of numerical methods, but in attuning to the fundamen- 
tal physics of the problem. The details of the algorithms 
are presented in the Appendix. 

II. QUANTUM-RING MODEL 

We focus on quantum rings realized in semiconduc- 
tor heterostructures, which can be modeled by localizing 
electrons to a 2D (xy) plane. We use the effective-mass 
approximation with GaAs parameters, i.e., the effective 
electron mass m* — 0.067 m e and the dielectric constant 
k = 12.7. The many-electron Hamiltonian is written in 
SI units as 

i JV n 2 

H = ^EHWi + eAfrif+V- £ - , 

2m* ^— i ^— ; Attcqk r, — r,- 

i=l i<j 

N 

+ ^2[V c An,0 t ) + V {mp {v l )+g*^ B Bs z ^. (1) 

The magnetic field B = Bz is chosen perpendicular to 
the xy plane, the vector potential is then, in linear gauge, 
A = —Bye x . The last term is the Zeeman energy that 
couples the external magnetic field with the electron spin. 
Here g* = —0.44 is the effective gyromagnetic ratio, 
Hb = eh/2m e is the Bohr magneton, and s z = ±^ for the 
up and down spins, respectively. The spin-orbit interac- 
tion is expected to be negligible in the GaAs structure 
having a wide band gap, it is therefore ignored in the 
Hamiltonian. The external potential that confines the 
electrons and defines the geometry of the quantum ring 
is chosen, in polar coordinates, to be 

V cxt (r,0) = \m*ulr 2 [1 + acos( P 9)] + V^l d \ (2) 

where hu>o = 5 meV is the confinement strength, and 
the cosine term defines the confinement geometry. We 
apply a circular and square shape determined by a — 
and (a,p) — (0.2,4), respectively^ The Gaussian term 
in Eq. |21 defines an antidot at the center, thus producing 
a ring-like shape of our system. We set Vb = 200 meV 
and the width parameter d = 10 nm, which is a sufficient 
value for two-dimensionality. The shapes of the total 
external potentials for the both geometries are shown in 

Fig.n _ 

Within the circular ring geometry (a = 0), we also ap- 
ply an impurity potential describing repulsive Coulom- 
bic impurities located randomly in the vicinity of the 
quantum-ring system. It is written as 

M = E 7 u \ W=V ' (3) 




FIG. 1: Shapes of the external confinement potentials deter- 
mining circular and square quantum rings with a = and 0.2 
in Eq.|5| respectively. 



where Mmp is the number of impurities, and Rj, and 
dk are their random lateral and vertical positions in the 
ranges of < Rk < 100 nm and < dk < 10 nm, respec- 
tively. This model, motivated by single-electron tunnel- 
ing experiments;^ has been applied recently in statistical 
studies on quantum dotsii^*i& Similarly to those studies, 
for each N[ mp = 5 ... 30 we apply 1000 spatial configura- 
tions in order to obtain good statistics. 



III. SINGLE-ELECTRON SPECTRA 

In order to obtain insight into the energy-level struc- 
ture in the quantum rings studied, we first computed the 
single-electron spectra of a system of noninteracting elec- 
trons as a function of the magnetic-field strength. Fig- 
ure El shows the eigenenergies £j for a circular (a) and 
square (b) quantum ring with up to i ~ 25 and B = 10 
T. The spectrum of our circular ring is similar to the one 
obtained using a pure parabolic 2D ring models Due to 
the finite ring width, the second Fock-Darwin level and 
the beginning of the third one can be seen in the spec- 
trum. Increasing the width by making d smaller brings 
the upper Landau levels lower in energy. 

In the spectrum of a square ring shown in Fig. |2{b), 
the four lowest energy levels are decoupled from the up- 
per levels and form a braid-like structure as a function of 
the magnetic field. Similar decoupling can be observed 
for the next four levels at low fields. This behavior re- 
sults from the four- fold symmetry of the square confining 
potential. The probability densities of the lowest eigen- 
states show that the four lowest levels correspond to ener- 
getically stable corner modes. The next four levels in the 
lower row, instead, correspond to the side modes. Valm- 
Rodrfquez and co-workers^ have analyzed modes of this 
type in their far-infrared-absorption studies on triangular 
and square quantum dots. 

In addition to the decoupled single-electron levels, 
there are several other avoided level crossings in the spec- 
trum of a square ring, particularly at low magnetic fields. 
The level repulsion in quantum dots is usually interpreted 
as a signature of quantum chaos£ Contrary to the square 
quantum dotji£ our square-ring system is non-integrable 
also at zero magnetic field. A similar system having steep 
walls (Sinai billiard) is a famous example of a chaotic 
system, it has been extensively studied in the context 
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FIG. 2: Single-electron spectra in circular (a) and square (b) 
quantum rings as a function of the magnetic field. The Zee- 
man splitting is omitted. 



1.5 



0.5 



E 









N imp = 5 
























■ 












- 




FIG. 3: Level-spacing distributions of noninteracting elec- 
trons at B — 10 T in impurity-affected quantum rings. The 
level spacings are calculated from one thousand random im- 
purity configurations for 5 (upper panel) and 30 impurities 
(lower panel). 



cases, leading to a two-peak structure. When the num- 
ber of impurities is increased to 30 (lower panel in Fig.|2J), 
the signatures of the shell structure clearly disappear. On 
a coarse scale, the distributions resemble Wigner-Dyson 
forms indicating a disordered system. However, the dis- 
tributions show fine structure resulting from the correla- 
tion between consecutive energy levels in the presence of 
a high number of impurities. 



of both classical and quantum billiardsiS In Fig. Efb) it 
can be seen that the number of avoided crossings de- 
creases at high magnetic fields because the electron or- 
bitals shrink due to the increasing magnetic confinement. 
However, the system remains chaotic, again in contrast 
with a square quantum dot which becomes integrable at 
B -> oo. 

The presence of external impurities (see Eq. |3J) lead 
to irregular deviations in the single-electron energy lev- 
els. In Fig. we plot distributions (gray scale) of single- 
electron level spacings at B = 10 T calculated from one 
thousand random impurity configurations for Aq mp = 5 
and 30, respectively. The level spacings correspond to 
the addition energies of noninteracting electrons defined 
as Ao(iV) = £at/2+i — cat/2- In the case of five impurities 
the distributions are approximately peaked around the 
level spacings of the clean quantum ring (cf. Fig. [5] at 
B = 10 T). However, when N — 14 (eg — £7), for exam- 
ple, there are two maxima in the distribution as seen in 
the upper panel of Fig- EI These result from the fact that 
while these two states are nearly degenerate at that field, 
they have a different slope with respect to B. Hence, 
in strongly disordered configurations the level splitting 
is considerably larger than in the approximately clean 



IV. MANY-ELECTRON PROPERTIES 

A. Circular and square geometries 

For the corresponding many-electron problem we apply 
the SDFT. Using the numerical scheme outlined in the 
Appendix, we have calculated the total energies of differ- 
ent spin states for TV = 1 ... 17 and for each magnetic- 
field strength up to 10 T in steps of 0.2 T. The ground 
state for each N is then defined as the spin state having 
the lowest energy E(N). Figure 01 shows the chemical 
potentials (jl(N) = E(N) - E(N - 1) for N = 10 ... 17 
in circular (upper blue lines) and square (lower red lines) 
quantum rings. We omit the low electron numbers in 
Fig. 01 to display the differences between two rings more 
clearly. The transitions in the ground-state spins are 
marked in Fig. 01 such that the filled circles denote an 
increase S — » S + 1 , and the open circles mark a decrease 
S — > S — 1. The crosses mark an increase of S — » S + 2. 
The rightmost points correspond to full spin polarization. 
The chemical potentials for N — 10, 11 and 12 in circu- 
lar ring agree with those of Emperador and co-workers. 18 
The trend of pairing between peaks for consecutive val- 
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FIG. 4: (color online). Chemical potentials for N = 10 ... 17 
in circular (upper blue lines) and square (lower red lines) 
quantum rings. The filled circles, open circles, and crosses 
correspond to the spin changes S — > S + 1, S — > 5 — 1, and 
S — > S + 2, respectively. 



ues of TV, as well as quasiperiodic oscillations analyzed 
in Ref. fl8l are found to continue up to higher electron 
numbers. We see the predicted violations in the pairing, 
e.g., at B ~ 1 and 2.2 T when N = 11 - 12. These 
are explained by Hund's rule, causing partial spin po- 
larization near the level crossings in the corresponding 
single-electron spectrum (see Fig. 01. 

As seen in Fig. 0J the evolution of /i(iV) as a function 
of B is qualitatively similar in circular and square rings, 
although the oscillations at low fields are considerably 
stronger in the circular case. This difference is due to 
the level repulsion of the square-shaped ring that leads 
to smoother behavior of the total energies. The tendency 
of the symmetry-breaking to even out the oscillations has 
been detected also in quantum dots containing external 
impurities The square ring shows also increased sta- 
bility as a small increase in n(N) when N = 4 or 8 [not 
plotted in Fig. [J]. This effect is presented and discussed 
in detail in connection with the addition energies below. 

The chemical potentials could be directly compared 
to Coulomb blockade oscillations of transport measure- 
ments, but we are not aware of such experiments for 
quantum rings. In the quantum dots, however, a remark- 
ably good agreement has been obtained between the ex- 
perimental result of the Coulomb blockade oscillations 9 
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FIG. 5: (color online). Addition energies (yellow: low, black: 
high) as a function of the magnetic held and up to 16 electrons 
in circular (a) and square (b) quantum rings. The open and 
filled circles mark the first and final full spin polarizations, 
respectively. The dashed lines are plotted to guide the eye. 



and SDFT calculations.^ 

Next we consider the second energy differences, i.e., the 
addition energies defined as A(./V) = n(N + 1) — n(N) = 
E(N - 1) - 2E(N) + E(N + 1). Besides chemical po- 
tentials, they are also measurable quantities^* and give a 
more detailed view on the energetic structure of quantum 
rings. Figure [3] shows the full phase diagram of A (AT) as 
a function of B and N for a circular (a) and square (b) 
ring. The light and dark regions correspond to low and 
high addition energies, respectively. The open and filled 
circles mark the first and final full spin polarizations, and 
the dashed lines are to guide the eye. The characteristic 
energy oscillations are evident in both geometries. The 
structure of the oscillations, however, shows interesting 
differences. In the circular case there are three easily dis- 
tinguished phases: (i) strong oscillations at low B that 
mostly correspond to spin changes — > 1 — > following 
Hund's rule; (ii) spin flip region at intermediate B that 
becomes broader as N increases; and (hi) a polarized 
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FIG. 6: Magnetization of a 24-electron circular quantum ring. 
Different magnetic phases characteristic of finite Fermion sys- 
tems are marked in the figure. 



regime in high magnetic fields on the right side of the 
dashed line(s), characterized by large- width oscillations 
of a period A<& ~ $ (the flux $ = BA, where A is the 
estimated ring area). A similar division can be found in 
the addition energies of a square ring [Fig. b)] . How- 
ever, the effect of the four-fold asymmetry is clear in the 
pronounced stability of the four-electron ring, ft results 
from the decoupling of the lowest single-electron levels 
shown in Fig. [5] The stability is strongest after the polar- 
ization at B ~ 3 . . . 4 T, since then all the four separated 
levels are (singly) occupied. Likewise, the eight-electron 
ring is particularly stable at low fields when S = 0. 

The addition energies in high magnetic fields are very 
similar between the circular and square rings, particu- 
larly when N is large. This result agrees qualitatively 
with the experimental observation of Fuhrer and co- 
workers^ It can be explained by the fact that the mag- 
netic confinement has a parabolic shape. Furthermore, 
the electron-electron interactions make the square ring ef- 
fectively more symmetric^ In the onset of the full spin 
polarization there are two differences. First, the square 
ring has a kink at N = 4 due to its specific energy spec- 
trum discussed in Sec. IIIII Secondly, there are many 
depolarizations in the square ring when N — 6 . . . 13, 
occurring between the first and second full spin polar- 
izations marked as open and solid circles in Fig. \E\ re- 
spectively. However, this effect is relatively faint and its 
emergence is very sensitive to the ring parameters. For 
example, depolarization can be obtained in the circular 
ring at N = 8 ... 11 if the Zeeman energy in Eq. 2] is 
slightly reduced. 

The different magnetic phases discussed above can be 
characterized in detail from the magnetization curves 
for large electron numbers. The magnetization is de- 
fined as the derivative of the free energy with respect 
to the magnetic field and reduces in zero temperature to 
M = —dEtot/dB. In Fig. we plot the magnetization 
of a 24-electron circular quantum ring as a function of 
the magnetic field strength B. In the two steps marked 



in the figure as (i) and (ii), the ground-state spin is 
zero and the (doubly) filled angular-momentum states are 
(l z ) = 0, —1, . . . , —11 and —1, —2, . . . , —12, respectively. 
They are calculated as the expectation values of the angu- 
lar momentum operator l z = —ih[x(d / dy) — y(d/dx)] for 
different KS states. The first occupation (i) from (l z ) = 
to —(N/2 — 1) directly corresponds to a quantum Hall 
state 20 with a filling factor // = 2. In the spin-flip region 
at B ~ 5 . . . 10 T the spin polarization of the ring grad- 
ually increases and the magnetization is characterized 
by short-period oscillations. In the first fully polarized 
(S = 12) state the occupation is (l z ) = —2, —3, . . . , —25. 
After that the oscillations become regular such that each 
step corresponds to an increase of 24 in the total angular 
momentum \L Z \ = | as the electrons jump from 

to (l^) — 1. Consequently, the hole in the electron 
density at the center of the ring increases, and there are 
no signs of edge reconstruction, which confirms the result 
of Emperador and co-workers 

The behavior in the electron occupations and mag- 
netization as a function of B is similar to that of a 
quantum-dot system^ The most obvious difference is 
that the v = 1 state with an occupation from (l z ) =0 to 
— (N — 1) is not found in quantum rings. This state, i.e., 
the maximum-density droplet^ is particularly stable in 
quantum dots and has been observed experimentally. 9 
However, we expect that in quantum rings the signatures 
of the v — 2 state and spin flips, as well as the regular 
oscillations in the polarized [y < 1) regime could be ob- 
served in magnetization experiments using, for example, 
sensitive micromechanical magnetometers^ 



B. Impurity effects 

Finally we analyze the effect of impurities on the ad- 
dition energies in a many-electron system. The applied 
impurity model is defined in Eq. [31 For simplicity, we 
focus here on the fully spin-polarized regime where the 
addition-energy oscillations in corresponding impurity- 
free quantum rings (see Fig. |3J) are regular. Figure [7| 
shows the addition-energy distributions of a 12-electron 
quantum ring at B = 10 T. For each number of impuri- 
ties we have calculated the electronic structure and ener- 
getics for one thousand random impurity configurations. 
We find that the distribution for Ni mp = 5 shows a clear 
two-peak structure. This is due to the fact that in this 
case the rings can be roughly divided into relatively dis- 
ordered and clean ones depending on the actual location 
of the impurities. These two types of rings have then 
different radii on the average, which eventually leads to 
a phase shift in the AB oscillations as a function of the 
magnetic field. 

As the number of impurities is increased, the two peaks 
gradually merge and finally at -/Vj mp = 30, correspond- 
ing to an impurity density of ~ 10 nm , we find a 
Gaussian-like symmetric distribution shown in Fig. EKd) . 
The symmetric shape is a result of the electron-electron 
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FIG. 7: Addition-energy distributions of a 12-electron quan- 
tum ring having 5 (a), 10 (b), 20 (c), and 30 (d) impurities. 
The magnetic field is set to B — 10 T corresponding to the 
fully-polarized regime. 



interactions that even out the tilted shape and irreg- 
ularities found in the level-spacing distributions (non- 
interacting electrons) shown in Fig. |21 The Gaussian 
shape is qualitatively similar to what has been found 
in disordered quantum dots both experimentallji 2 ^ 2 ^ and 
theoretically using Hartree-Focl*2& and density- functional 
methodsii§i2L2£ In addition, modifications to the ran- 
dom matrix theory that take the interactions into ac- 
count have led to a good agreement between the exper- 
iments and the theory^ whereas the bare random ma- 
trix theory yields the Wigner-dyson distribution. Hence, 
we expect that the distribution shown in Fig. |7fd) for 
quantum rings corresponds to the limit of disorder which 
ultimately indicates the disappearance of the AB oscilla- 
tions. Following the analogy to quantum-dot systems, 28 
we expect that a Gaussian-like distribution could be ob- 
tained in conductance experiments for large (N > 100) 
quantum rings. In that case the addition-energy statis- 
tics would be determined as a function of N. 



V. SUMMARY 

We have used efficient computational algorithms 
for solving the Kohn-Sham equations of spin-density- 
functional theory to examine the effects of geomet- 
ric deviations on the ground-state properties of two- 
dimensional quantum rings in magnetic fields. Due to the 
circular shape of the magnetic confinement, increasing 
magnetic field evens out the variations between the ener- 
getics of quantum rings of different geometries. The mea- 
surable quantities in both systems can be characterized 



by quasiperiodic oscillations, of which we have been able 
to identify the quantum-ring counterparts of the integer 
quantum-Hall states. However, the Aharonov-Bohm os- 
cillations at high fields are very sensitive to the presence 
of external impurities which may produce a systematic 
phase shift. In the high- impurity (disorder) limit the 
electron-electron interactions make the addition-energy 
distributions Gaussian-like. This behavior, indicating 
the disappearance of the oscillations, is similar to what 
has been found in quantum dots. 



APPENDIX: NUMERICAL SCHEME 

Our numerical procedure for solving the electronic 
properties of the system within the SDFT consists of two 
parts, namely, the solution of an effective Schrodinger 
equation [the Kohn-Sham (KS) equation] and the self- 
consistent determination of the spin-densities. The KS 
equation in a magnetic field is 



1 



2m 



-{-ihV + eA{r)Y + V KS (r) 



*f (r) = e^Hr), 



(A.l) 

where the KS potential Vks(r) is a sum of the exter- 
nal potential defined above, the Hartree potential, and 
the exchange-correlation potential given as V^ c (r) = 
6E xc [p, £]/5p' T (r). Here p a are the electron spin densities, 
tr denotes the spin index, and £(r) = [p^ (r) — p J '(r)]/p(r) 
is the local spin polarization. For E xc we use the local 
spin-density approximation with the functional provided 
by Attaccalite and co-workersiSi 

The lowest n solutions of the eigenvalue problem l|A.l|) 
are obtained by applying the evolution operator, 



-tH 



(A.2) 



repeatedly to a set of states {ipj, 1 < j < n}, and or- 
thogonalizing the states after every step. Instead of the 
commonly used second-order factorization in combina- 
tion with the Gram-Schmidt orthogonalization, we use 
the fourth-order factorization for the evolution operator 
given by^i^ 



where 



T —-tV - 



(A.3) 



T 



1 



2m* 



-II 2 



2m* 



n 2 ,], n = p + eA(r) 



(A.4) 

is the kinetic-energy operator. We have defined the local 
modified potential 3 "^^ as 



V = V+-e 2 [V,[T,V]]. 



(A.5) 



Note that the vector potential A(r) does not contribute 
to the commutator. 
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We have shown in Ref.^]that the density matrix e~ eT 
can be exactly decomposed for a uniform magnetic field 
as 



where £ = eheB/m*, and 



(A.6) 



COSh(e) - 1 and C M (£) = ^. (A.7) 



£sinh(0 



The exact factorization shown in Eq. (|A.6() is possible 
because the Hamiltonian of an electron in a uniform mag- 
netic field is quadratic, and higher-order commutators of 
TL 2 and Tly are either zero or simply proportional to H 2 
and Hy. The two key commutators are 



[n?,prLn?]] = 8h 2 e 2 B 2 n 2 l . 



(A. 



Hence, all higher order commutators appearing in the 
Baker-Campbell-Hausdorff formula can be summed back 
to the original operators H 2 and H 2 . 

The above result can be easily generalized to a charged 
particle in a uniform magnetic field in an arbitrary ex- 
ternal potential: Inserting the exact factorization 1A.6|) 
into the factorization (|A.3(I of the full Hamiltonian yields 
the final result: 



x e 4 " 



x e 4 " 



rC B (f)n= e -f £ v e - 



(A.9) 



rC M (|)n= 



e 



0(e 5 



we have shown that this algorithm is, depend- 
ing on the system and the desired accuracy, a factor of 
10 to 100 more efficient than the second order factoriza- 
tion. There, we have worked in circular gauge, here we 
point out that the method becomes even more efficient 
in linear gauge. The algorithm is then applied as follows: 

(1) Start with a set of suitably-chosen initial states in 
the coordinate space. 



(2) Multiply these states by e 



(3) Fourier transform the x-coordinate of each state to 
the fc x -space and multiply by e _ 4^T c ' E (l)(' c ^+ eS a) 2 . 

(4) Fourier transform now y to the fc y -space, and mul- 
tiply by e~^ CM ^ k2 y. 

(5) Do the inverse transformation back to y and mul- 
tiply by e~4^ c E(l)( fe -+l B y) 2 . 

(6) Fourier transform k x back to the x-space and mul- 



tiply by e" 



Then repeat the steps (3)-(5) and 



finally multiply the states by e" 



(7) Orthonormalize the states and repeat the proce- 
dure until convergence has been obtained. 



Thus, the implementation of the algorithm requires the 
equivalent of two 2D Fourier transforms. In other words 
it is computationally no more costly than the case with- 
out a magnetic field. For the orthonormalization, we di- 
agonalize the matrix of the overlap integrals and from 
these we construct a new set of orthonormal states. 

A persistent problem of density-functional calculations 
is that the naive charge density mixing iteration scheme, 
usually require a large number of iterations for conver- 
gence. We overcome this problem by using a methodM 
which solves for p(r) directly by applying a Newton- 
Raphson procedure is used. We define 



Ap CT (r) 



h 



n CT (h)|$n PT , Pi ](r)| 2 -p CT (r) (A.10) 



as the density difference between two self-consistent it- 
erations. Here n 17 is the occupation factor, are the 
orthogonalized solutions of Eq. i|A.l|) . and p a (v) is the 
density used for the calculation of Vks- The sum goes 
over all occupied (hole: h) states. Then the density cor- 
rection Sp a (r) is determined by a linear equation 



dV e a ' a ' 



(r,r';0)£p ff (r'). (A.ll) 



Here d is the dimension of the system and e a ' a is the 
static dielectric function of a non-uniform electron gas^ 
It contains the zero-frequency Lindhard function and the 
particle-hole potential V^{v,v') = 5V£ s (r) / 6 p°' (r'). 
To avoid the calculation of unoccupied (particle: p) 
states we seek for an approximation for the static re- 
sponse function that only needs the calculation of oc- 
cupied states. For the purpose such an algorithm, we 
recall that linear response theory can be derived^ from 
an action principle for excitations of the form \ij)(t)) = 
ex P(J2ph c ph(t)ala h ) |0o)j where \<f> ) is the ground state, 
and c p h(t) are particle hole amplitudes. If we assume 
that the particle-hole amplitudes are matrix elements of 
a local function uo a (r), i.e., c p h{t) = (p\ w CT (r, t) \h) we end 
up with Feynman's theory of collective excitations^ Us- 
ing the commutator of the kinetic part of the effective 
Schrodinger equation with u a , 



[n 2 ,< 



±[n [n,^] + [n,^]n] 
-|[n(Vw ff ) + (v^)n], (A.12) 



the magnetic part completely cancels out because it is lo- 
cal. Thus, applying the response- algorithm in a magnetic 
field also applies no computational overhead compared to 
the zero-field case. In the "collective approximation" we 
can rewrite Eq. IA.11I asi£ 



±V • [/(r)V] + 2£ CT , S° F * V%?1 * Si * 

where now 

&p°(v) = A/(r)-^(r,rV» ff (r'), 



(A.13) 
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and 
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